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Abstract. We prove that the Leech lattice is the unique densest lattice in R'^*. 
The proof combines human reasoning with computer verification of the prop- 
erties of certain explicit polynomials. We furthermore prove that no sphere 
packing in R-^* can exceed the Leech lattice's density by a factor of more than 
1 + 1.65 ■ 10"^'^, and we give a new proof that Eg is the unique densest lattice 
in R**. 



It is a long-standing open problem in geometry and number theory to find the 
densest lattice in R". Recall that a lattice A C M" is a discrete subgroup of rank n; 
a minimal vector in A is a nonzero vector of minimal length. Let |A| ~ vol(R"/A) 
denote the covolume of A, i.e., the volume of a fundamental parallelotope or the 
absolute value of the determinant of a basis of A. If r is the minimal vector length 
of A, then spheres of radius r/2 centered at the points of A do not overlap except 
tangentially. This construction yields a sphere packing of density 



since the volume of a unit ball in R" is 7r"/^/(n/2)!, where for odd n we define 
{n/2)\ = r(n/2 + 1). The densest lattice in R" is the lattice for which this quantity 
is maximized. 

There might be several distinct densest lattices in the same dimension. For 
example, the greatest density known in R^^ is achieved by at least 23 distinct 
lattices, although they are not known to be optimal. (See pages xix and 178 of 
[CS2j for the details.) Wc will speak of "the densest lattice" because it sounds 
more natural. 

The problem of finding the densest lattice is a special case of the sphere packing 
problem, but there is no reason to believe that the densest sphere packing should 
come from a lattice. In particular, in R^° the densest packing known is the Best 
packing Pioc, which is not a lattice packing (see [CS2[ p. 140]), and it is conjectured 
that lattices are suboptimal in all sufficiently high dimensions. However, many of 
the most interesting packings in low dimensions are lattice packings, and lattices 
have strong connections with other fields such as number theory. For example, 
the Hermite constant 7„ is defined to be the smallest constant such that for every 
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positive-definite quadratic form Q{xi, . . . , Xn) of determinant D, there is a nonzero 
vector (ui, . . . , u„) G Z" sucli that Q{vi, . . . , -Wn) < JnD^^'\ Finding the maximum 
density of a lattice packing in M" is equivalent to computing 7„. 

The densest lattice in M" is known for n < 8: the answers arc the root lattices 
Ai, A2, A3, 1)4, 1)5, Eq, E-j, and E^. For n = 3 this is due to Gauss [G], for 
4 < n < 5 to Korkine and ZolotarefF [KZn[KZ2] . and for 6 < n < 8 to Blichfeldt 
[BI] . However, before the present paper no further cases had been solved since 1935. 
In 1946 Chaundy claimed to have dealt with n = 9 and n ^ 10, but his paper [Ch| 
implicitly assumes that a densest lattice in K" must contain one in R"~^ as a cross 
section. That is known to be false (see [CSlj l. so the paper appears irreparably 
flawed. 

In each of the solved cases, the optimal lattice is furthermore known to be unique, 
up to scaling and isometrics. This was proved simultaneously with the optimality 
for ri < 5, for n = 6 it was proved by Barnes [Ba| . and for 6 < n < 8 it was proved 
by Vetcinkin JVeJ. 




In this paper we deal with n = 24 (the theorem numbering is as it will appear 
later in the paper): 

Theorem 19. 3i The Leech lattice is the unique densest lattice in R^"*, up to scaling 
and isometrics ofM.'^^. 

In terms of the Hermite constant, 724 — 4. We also give a new proof for Es'- 

Theorem 111.7] (Blichfeldt. Vetcinkin). The Eg, root lattice is the unique densest 
lattice in R^, up to scaling and isometrics ofM.^. 

Our work is motivated by the paper |CE| by Cohn and Elkies (see also |Co) ). 
which proves upper bounds for the sphere packing density. In particular, the main 
theorem in [CEj is an analogue for sphere packing of the linear programming bounds 
for error-correcting codes: given a function satisfying certain linear inequalities one 
can deduce a density bound. It is not known how to choose the function to optimize 
the bound, but in and R^'' one can come exceedingly close to the densities of 
Es and the Leech lattice, respectively. This observation, together with analogies 
with error-correcting codes and spherical codes, led Cohn and Elkies to conjecture 
that their bound is sharp in R* and R^'*, which would solve the sphere packing 
problem in those dimensions. While we cannot yet fully carry out that program, 
in this paper we show how to combine the methods of [CEj with results on lattices 
and combinatorics to deal with the special case of lattice packings. We will deal 
primarily with the Leech lattice, because that case is new and more difficult, but 
in Section [n] we will discuss £^8- 

One might hope to use a relatively simple method. Section 8 of |CE| shows how 
to prove that the Leech lattice is the unique densest periodic packing (i.e., union 
of finitely many translates of a lattice) in R^*^, if a function from R^^ to R with 
certain properties exists. Using a computer, one can find functions that very nearly 
have those properties, and the techniques from Section 8 of [CEj can then be used 
straightforwardly to prove an approximate version of the uniqueness result: every 
lattice that is at least as dense as the Leech lattice must be close to it. It is known 
(see |Ma| p. 176]) that the Leech lattice is a strict local optimum for density, among 
lattices. Thus, if one can show that every denser lattice is sufficiently close, then it 
proves that the Leech lattice is the unique densest lattice in R^**. 
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Unfortunately, this approach seems completely infeasible if carried out in the 
most straightforward way. When one naively imitates the techniques from Section 8 
of [CEj in an approximate setting, one loses a tremendous factor in the bounds, 
and that puts the required computer searches far beyond what we are capable of. 
In this paper we salvage the approach by using more sophisticated arguments that 
take advantage of special properties of the Leech lattice. In particular, we make 
use of three beautiful facts about the Leech lattice: its automorphism group acts 
transitively on pairs of minimal vectors with the same inner product, its minimal 
vectors form an association scheme when pairs are grouped according to their inner 
products, and its minimal vectors form a spherical 11-design. (Note that the second 
property follows from the first, but our work uses another proof of it, from [DGSj .) 

Our proof depends on computer calculations in some places, but they can be 
carried out relatively quickly, in less than one hour using a personal computer. The 
calculations are all done using exact arithmetic and are thus rigorous. We have 
fully documented all of our calculations and made available commented code for 
use in checking the results or carrying out further investigations. See Appendix \X\ 
for details. 

Appendix IB] contains very brief introductions to several topics: the Leech lattice, 
linear programming bounds, spherical designs, and association schemes. For more 
details, see |CS2] . The expository article [E2J also provides useful background and 
context, although it does not include everything we need. 



2. Outline of proof 

We wish to show that the Leech lattice, henceforth denoted by A24, is the unique 
densest sphere packing among all lattices in M^^. Let A be any lattice in M^^ that is 
at least as dense as A24. Without loss of generality, we assume that A has covolume 
1. Then the restriction on its density simply means its minimal vectors have length 
at least 2. 

We first show, using linear programming bounds, that A has exactly 196560 
vectors of length approximately 2 (called nearly minimal vectors), and that the 
next smallest vector length is approximately y/6. 

We rescale those 196560 nearly minimal vectors to lie on the unit sphere. Then 
they form a spherical code with minimal angle at least (p, where cos f is very near to 
(and greater than or equal to) 1/2. Note that in S^^ there is a unique spherical code 
of this size with minimal angle Tr/3 — cos~^(l/2), and it is the kissing arrangement 
of A24; the spherical code derived from A should be a small perturbation of this 
configuration. 

Using linear programming bounds, we show that the inner products of the unit 
vectors are approximately 0, ±1/4, ±1/2, ±1. We prove that if pairs of vectors are 
grouped according to their inner products, then they form an association scheme 
with the same valencies and intersection numbers as in the case of A24, and that 
it must therefore be the same association scheme. This isomorphism gives us a 
correspondence between minimal vectors of A24 and nearly minimal vectors of A, 
such that corresponding inner products are approximately equal. 

Using this correspondence, we find a basis of A whose Gram matrix is close to 
the Gram matrix of a basis of A24. Finally, from the strict local optimality of the 
Leech lattice we conclude that A must in fact be the Leech lattice. 
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3. Notation 

We begin by recording our normalizations of some special functions (which are 
always as in [A AR] ) . and defining some notation. 

The Laguerre polynomials Lf{z) are defined by the initial conditions Lq{z) = 1 
and Li(z) = 1 + a — z and the recurrence 

iLf{z) = {2i-l + a- z)Lf_^{z) - {i + a ~ 1)L?_2(^) 

for i > 2. They are orthogonal polynomials with respect to the measure e~^x°' dx 
on [0, oo). If a = n/2— 1, then the functions on M" given by x q-'^M Lf {2'k\x\'^) 
form an orthogonal basis of the radial functions in L^(M"), and they are also eigen- 
functions of the Fourier transform with eigenvalue (—1)* (see (4.20.3) in [Lebj ) . 
Here we normalize the Fourier transform by 

m = / /(x)e2-<--*) dx. 

Note also that with this normalization of the Fourier transform, the Poisson sum- 
mation formula states that 

ksA ' ' teA* 

if / : M" ^ M is a Schwartz function, A C K" is a lattice, and 

A* = {y e R" : (x, y) e Z for all x e A} 

is its dual lattice. (See (28) in [H p. 49].) 

The ultraspherical (or Gegenbauer) polynomials C^{z) are defined by the initial 
conditions Cq (z) = 1 and Ci{z) ~ 2Xz and the recurrence 

zC^iz) = 2(* + A - l)zC^i(z) - (z + 2A - 2)C^_^iz) 

for i > 2. They are orthogonal polynomials with respect to the measure 

(1_:j2)A-1/2 

on [—1,1]. When A = n/2 — 1, that measure is proportional to the projection of the 
surface measure from onto an axis, and the ultraspherical polynomials play 

a fundamental role in the theory of spherical harmonics in M". Up to scaling, the 
ultraspherical polynomial C^' is the same as the Jacobi polynomial Pj;°''"\ where 
a = A- 1/2. 

Whenever we use Laguerre or ultraspherical polynomials, we will always set 
a = A = n/2— 1, where n = 24 in the Leech lattice proof and n = 8 in the Es 
proof. The term "ultraspherical coefficient" will mean a coefficient in the expansion 
of a polynomial as a linear combination of ultraspherical polynomials. 

Throughout this paper, A24 will denote the Leech lattice, and A will denote any 
lattice in M^'' that is at least as dense and satisfies |A| = 1 (except in Section [51 
where A denotes an arbitrary lattice, and in Section [11] which deals with Eg)- We 
think of A as being an optimal lattice, but we will not use that assumption. We 
do not even need to know a priori that a global optimum for density is achieved, 
although [CjL[ §17.5] shows that it is. 

Whenever /: R" ^ R is a radial function and r £ [0, 00), we will write /(r) for 
the common value f{x) with x E R" satisfying |a;| — r. 
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The surface volume of the unit sphere S*" ^ C M" will be denoted by 

It is important to keep in mind that vol(5"^^) is not the volume of the enclosed 
ball. 

We will make use of the numerical values 

£ = 6.733 • 10 
/i = 3.981 • 10 
V = 3.219 • 10 
uj = 1.703 • 10 

throughout the paper. Each will be defined the first time it is used, but we have col- 
lected the values here for easy reference. Note that terminating decimal expansions 
such as these represent exact rational numbers, not floating point approximations. 

4. Nearly minimal vectors 

In this section, we show that A must have exactly 196560 vectors of length near 2 
(this will be made more precise below) . The first subsection examines which vector 
lengths are possible in A, and the second then counts the nearly minimal vectors. 

4.1. Restrictions on the lengths of vectors. Let /: R be a radial 

Schwartz function with the following properties: /(O) — /(O) = 1, f{x) < for 
|a;| > r (for some number r), and /(t) > for all t. Proposition 3.2 of [CE| says 
that if such a function exists, then the sphere packing density in is bounded 
above by 

m V2/ ■ 

If we could find such a function with r — 2, then it would prove that A24 has the 
greatest density among all sphere packings in M.'^^, not just lattice packings. Cohn 
and Elkies conjecture that such a function exists, but the best they achieve in [CE] 
is r < 2 • 1.00002946. 

Our proof begins by constructing an explicit function / with 

r< 2(1 + 6.851- 10^32). 

Note that the existence of such a function proves that no sphere packing in can 
have density greater than 1 + 1.65 • lO"'^" times the density of A24. 

Unfortunately, the function we construct is extremely complicated (it would take 
far too much space to write it down here). It consists of a polynomial of degree 
803 with rational coefficients, evaluated at 27r|xp and multiplied by e~'^'^' . It 
was constructed by a lengthy computer calculation to optimize the value of r using 
Newton's method, and even verifying that it has the properties used below requires 
a computer, although fortunately that is much easier than finding the function. The 
accompanying computer file verifyf .txt includes code to verify all the assertions 
about / in this subsection of the paper. See Appendix [K\ for details. 

We can use techniques similar to those in [CE| to study the size of the short 
vectors in A. First, we need two lemmas. Set 

£ = 6.733 • 10^2^ 



-27 
-13 

and 

-11 
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and call a nonzero vector in A nearly minimal if it has length at most 2(1 + e). The 
reason for this choice of e will be apparent from Proposition 14. 31 below. 

Consider what happens if we rescale the nearly minimal vectors so that they all 
lie on the unit sphere. These vectors determine a spherical code 

Ca = {u/\u\ : u a nearly minimal vector} 

on the unit sphere S"^^, and the following lemma bounds its minimal angle. (See 
Appendix IB] for background on spherical codes.) 

Lemma 4.1. If u and v are nearly minimal vectors with u v, then the angle ip 
between u and v satisfies 

^°^^^^-2(TT7F- 

Proof. We have G [2. 2(1 + e)] and |m — wl > 2. By the law of cosines, 

\u\^ + \v\^ ~ \U ~ V\^ |u|2+|w|2_4 
COS ip = -; — —. < — — —. . 

2|u||i;| - 2|u||w| 

The bound (|up + |fp— 4)/(2|u||'i;|) is convex as a function of | u | and | v \ individually, 
and hence it is maximized at one of the vertices of the square [2, 2(1 + e)]'^. In fact, 
the maximum occurs when \u\ = \v\ = 2(1 + e), in which case the bound becomes 

1 



1 - 



2(1 +6)2 • 

□ 



Lemma 4.2. There are at most 196560 nearly minimal vectors in A. 

Proof. This lemma is a straightforward application of the linear programming 
bounds for spherical codes (see Chapter 9 of [CS2] , or Appendix [B] for a brief 
summary). Let 

/e(.) = i^e(x + 1) (x + i) ' (x + i) ' .2 _ 1 ) ^ _ _ 

where the constant is chosen so that has zeroth ultraspherical coefficient 1. 
(The normalization is irrelevant for this proof, but it will be important later in the 
paper, so we use it here for consistency.) If s were 0, this polynomial would be the 
one used to solve the kissing problem exactly in (ggg Chapter 13 of |CS2] ). With 
the current value of e, the polynomial f^ has nonnegative ultraspherical coefficients 
and still proves that there are fewer than 196561 spheres in any spherical code in 
IR24 with minimal angle as in Lemma ITTl (We check this assertion in the computer 
file verifyrest.txt. In fact, the bound is less than 196560+ 10~^^.) Thus, there 
can be at most 196560 nearly minimal vectors in A. □ 

In addition to the definition of e above, set 

fi = 3.981 • 10"^^ 

1/ = 3.219 • 10"^^ and 

uj = 1.703 • 10"". 

Proposition 4.3. Every nonzero vector in A has length in 

[2,2(l + e)) U (\/6(l-/i),V6(l + /i)) U (V8(l - z^), V8(l + i^)) U (710(1 - tj), oo) . 
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Proof. By the Poisson summation formula, 

xeA teA* 

Because A is at least as dense as A24 (and has covolume 1), all nonzero vectors 
X G A satisfy \x\ > 2. Because r < 2(1 + e), by Lemma [4.21 there can be at most 
196560 vectors in A with 2 < |a;| < r. Within that range, / is a decreasing function 
of the radius, and we have 

196560/(2) < 1.644104221 • 10"^° 

(recall from Section [3] that /(2) denotes the common value f{x) with \x\ = 2). 

The key properties of / are /(O) = /(O) = 1, f{x) < for > r, and f{t) > 
for all t. It follows that 

1 + 1.644104221 • 10"^° + ^ f{x) > 1, 

where Af is the set of vectors in A at which / is negative. No vector in A can occur 
within any region on which / is less than (l/2)(-l. 644104221) • 10^^"; the extra 
factor of 2 comes from the fact that f{—x) = f{x) since / is a radial function (if 
X £ J\f then — X G J\f as well). That rules out all radii in the set 

[2(1 + e), V6{1 - fi)] U [V6(l + /i), %/8(l - ly)] U [V8{1 + z^), VT0(1 - uj)] . 

In the computer file verifyf .txt we prove this by examining the radial derivative 
of/. □ 

Much of the rest of the proof would still work if e, /i, v, and lu were somewhat 
larger. The two main places where they must be small are the final inequality (|9.4|) 
and the intersection number calculations in Subsection 16.21 (as well as the bounds 
used there). In each case they could be slightly larger, but not by a factor of 100. 

4.2. 196560 nearly minimal vectors. We can now show that there are exactly 
196560 nearly minimal vectors. We know from Lemma [4.21 that there are at most 
196560 of them. For the other direction, a lower bound greater than 196559, we 
need a new kind of linear programming bound. Recall that we have shown that 
all nonzero vectors in A are either nearly minimal or have lengths greater than 
\/6(l-/i)- 

Suppose we knew that all nonzero vectors either have length exactly 2 or at least 
^/E. (We will first explain our method under these overly optimistic hypotheses. 
Lemma 14.41 will then apply it using the actual bounds we have proved.) One might 
hope to count the nearly minimal vectors using a Schwartz function g: R^** ~* M 
such that g{x) < for |a;| > \/6, g{t) > for all t, and g{2) > 0. Given such a 
function, Poisson summation implies that 

x&A teA* 

and hence 

g{0) + Ng{2) > g{0) 

if there are N minimal vectors. Thus, N > {g{0) — g{Q))/ g{2). We conjecture that 
g can be chosen so that (g(0) — g{0j)/g{2) = 196560, which is the largest possible 
value because the Leech lattice has 196560 minimal vectors. We can construct 
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functions that come quite close to this bound, and wiU use one of them to prove 
the following lemma. 

Lemma 4.4. There are more than 196559 nearly minimal vectors in A. 

Proof. Define zi,. . . , zio by Zi = [47r(i + 1)10^J/10^. In other words, they have the 
following values: 



i 


1 


2 


3 


4 


5 


Zi 


25.13274122 


37.69911184 


50.26548245 


62.83185307 


75.39822368 



i 


6 


7 


8 9 


10 


Zi 


87.96459430 


100.53096491 


113.09733552 125.66370614 


138.23007675 


There are unique 


coefficients ai, 


. . . , asr such that 








1- 


37 





i=l 

has a single root at Z2 and a double root at Zi for i > 3, and 

37 

l + J2i-iya,Ll\x) 

i=l 

has a double root at each Zi for i > 1. Neither polynomial has any other nonnegative 
roots. (We check this in the computer file verifyg.txt using Sturm's theorem, 
except that we do not check the uniqueness of the coefficients because we do not 
require it.) 

Define 5 : 1^24 ^ ]g 

It follows that 

= (l + f2i-iya,Ll\27r\x\')^ e-^^^^ 

We have g{x) < for \x\ > V6{1 — /i), because Z2 < 2tt ■ 6(1 — n)"^, the function g 
changes sign only at Z2, and g{0) > 0. For all t, we have g{t) > 0. 
Apply Poisson summation to g, to deduce 

xeA tsA* 
Applying the two inequalities above shows that 

g{0) + 9i^) > 9{0), 
xeM 

where A4 is the set of nearly minimal vectors. The function g{x) is positive and a 
decreasing function of \x\ on the interval [2, 2(1 + s)], so 

5(0) + \M\g{2) > 5(0). 

However, 

> 196559, 

5(2) 
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SO there are more than 196559 nearly minimal vectors. (All these inequalities are 
checkedinverifyg.txt.) □ 

Thus, there must be exactly 196560 nearly minimal vectors, as desired. We 
conjecture that this method could be used to recover each of the coefficients of the 
Leech lattice's theta series, but we will not need that for our proof. 



5. Inner products in the spherical code 



We now continue to study the polynomial 



1\ / 1\ o A 1 



from the previous section. 
Note that 



- 1 i 

2(l + e)2 < 2 



Thus, all the inner products except 1 in the spherical code Ca are at most 1/2 + e. 
We seek bounds for how far from 0, ±1/4, ±1/2, ±1 they can be. The ±1 cases 
must be exact, because of Lemma l4.ll and the fact that u G Ca iff —u G Ca (i.e., 
the code is antipodal). 

Because the zeroth ultraspherical coefficient of /e is 1, it follows from the usual 
proof of the linear programming bounds for spherical codes (see Appendix [B| that 

196560/,(1) + J2 feii^^ y)) > 196560^, 

where the sum is over vectors in the spherical code. Because the four inner products 
{x,y), {—x,—y), {y,x), and {—y,—x) are all equal, and all the terms in the sum 
are nonpositive, we see that no term in the sum can be less than 

196560^ - 196560/e(l) 
4 ' 

which is approximately —1.99 • 10~^^. 

Now a short calculation implies that all the inner products must be within 6.411 • 
10~^ of one of the numbers 0, ±1/4, ±1/2, ±1. The exponent is only 9 because has 
double roots, and one can stray quite far from a double root without substantially 
changing the function's value. 

In the rest of the paper, we make the following definition: let 

<j — max min{|(u, w) — a| : a G 5}, 

where S = {0, ±1/4, ±1/2, ±1}. In other words, a is the maximum "error" in the 
inner products. We have just shown that a < 6.411 • 10~^. We will improve the 
upper bound for a substantially, and ultimately we will show a — 0. 

5.1. Better bounds for a. Recall that we computed in Section WA\ that vectors 
of A with length close to must have length in the interval V6{1 — fj.), V6{1 + /i). 
Therefore if we have nearly minimal vectors u, v with {u, w) « 1 (with error less 
than 10~^, say), then we see that |u — wj ~ \/6. Therefore 

6(1 - fJ.f < {u, u) + {v, v) - 2{u, v) < 6(1 + fif. 
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In addition we have 

4 < {u,u) < 4(1 + e)^ 



Therefore 



so 



- 6(1 + fif < 2{u, v) < 8(1 + ef - 6(1 - fif, 
4-3(l + Ai)2 ^ / u v\^,, , .2 3, ,2 



4(1+.)^ -^^'-'^^ 

Similarly we get for (u, v) ~ that 

1- {l + uf I u V 



(l + £)2 
and for {u, v) k, 2 that 



2 - (1 ^ / u V \ ^ ^ 2 1 



because 2 < |u — < 2(1 + e). Combining these (and if necessary replacing v 
with — w), we find that for u ^ ±z;, the inner product differs from an 

element of {0, ±1/4, ±1/2} by at most 6.43801 • 10"^^ We conclude that the error 
a in the inner products is at most 6.43801 • 10^^^. 

However, we will be able to get much better bounds in Section [71 once we have 
shown that the spherical code Ca gives us an association scheme. 

6. Association schemes 

We would like to turn the 196560 points in the spherical code Ca into a 6-class 
association scheme by grouping pairs according to their approximate inner products. 
(See Appendix iBl for background on association schemes.) It is not clear that this in 
fact defines an association scheme, but we will show that it does. Furthermore, we 
will show that this association scheme is isomorphic to the one derived from A24. To 
achieve this, we show that the intersection numbers are the same as in A24. That 
will also show that it is an association scheme, by showing that the intersection 
numbers are independent of the pair of points. We use the same techniques as 
[DGS| . but we need to keep track of error bounds. 

6.1. Spherical design. First, we show that Ca is nearly a spherical 10-design. 
(See Appendix IB] for background on spherical designs.) Let 

Cll(l) ' VOl(523) ■ 

The advantage of this normalization of the ultraspherical polynomials is that for 
every finite subset C of 5'^^, 



E C'.((x,z/)) 



Eev,(z) 



zee 



where evi(z) denotes the evaluation at z map in the dual space to the i-th degree 
spherical harmonics. Although this fact is well known (e.g., it is equivalent to 
Theorem 9.6.3 in [A AR| ) . we will explain it here for completeness, because the 
correct normalization is important for our application. 
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Lemma 6.1. IfC is a finite subset of S^^, then 



x,y£C 



zee 



Proof. Let d be the dimension of the space of spherical harmonics of degree i, and 
let S\,...,Sd be an orthonormal basis of that space. By Theorem 9.6.3 of [AARj . 

d 

Y,SjHS,iz)^a{{w,z)). 
Let f = o-jSj be any spherical harmonic of degree i. Then 



(ev,(z))(/) = ^a,5,(z) 

d ^ 

= Y] Sj (z) / Sj {w)f{w) dw 



{llS,{w)S,{z)^ f{w)dw 



523 



Ci{{w, z))f{w) dw. 



Thus, 



C,{{w,z)) f{w)dw. 



In other words, applying the element X^^ec of the dual space is the same as 

taking the inner product with 



w 



EQ(Kz)). 



It follows that 



Eevi(z) 



zee 



= V / C^{{w,x))C^{{w,y))dw 



x,yee ' 

= ^{e^i{x)){w^Ci{{w,y))) 

x,yeC 
x,yeC 

as desired. 

The following lemma asserts that Ca is nearly a spherical 10-design: 



□ 
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Lemma 6.2. Ifg: S^^ 



zGCa 



M is a polynomial of total degree at most 10, then 
196560 



V0l(S'23) 



5.23 



g{z) dz 



< 2.50193- 10"^|g|2, 



where \g\2 denotes the norm on L [S ). 

Proof. Without loss of generality we may assume that g is a harmonic polynomial 
(for every polynomial on M^^, there is a harmonic polynomial of equal or lesser 
degree with the same restriction to S*^"^; see equation (5) in [HI p. 17]). We will use 
the polynomial 



f,{x) = K,{x + l)[x + 



X^ \ X 



1 



2(1 + e)^ 



from earlier in the paper. Recall that is normalized to have zeroth ultraspherical 
coefficient 1 (using the standard normalization of the ultraspherical polynomials, 
not the new normalization d). If Ci denotes the coefficient of d in /e, then 



196560/,(1)> fe{{x,v)) 



10 



1965602 + ^c, ai{x,y)) 

4—1 x,y£C 



10 



= 196560^ + ^0, 



i=l 



zGCa 



from which it follows that 



10 

1=1 



< 196560/e(l) - 196560^ < 7.9775 • 10" 



Therefore, 



10 

E 



E ^^^(^^ 



i=l 

We can now bound 



< (7.9775 • 10"^^) max - < 6.25964 • 10"^°. 

» C,; 



V- , ^ 196560 f 

1^ 5^ T7J23T / 



g(z) dz 



Write g = X]i=o52' where gi is homogeneous of degree i. The integral cancels with 
the 5o term in the sum, so we simply need to bound 

10 

E E 

i=l zGCa 
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For that, we use the definition of the norm and the Cauchy-Schwarz inequahty to 
deduce that 



10 

<E 



E ^^*(^) 



i=l z£Ca i= 


=1 zeCA 








10 




2 


10 






E ^^'(^) 


\ 


El5»l2 




i=i 


z&Ca 







< 2.50193- 10-^|5|2, 



as desired. 



□ 



In fact, it follows immediately that Ca is nearly a spherical 11-design in the same 
sense, because it is antipodal (if x e Ca then —x e Ca) and thus every homogeneous 
polynomial of odd degree averages to over Ca- However, we will not need that 
fact. 

It is worth pointing out for completeness that the minimal vectors in A24 do not 
form a spherical 12-design: if y is a minimal vector then the polynomial 

(16 - {x, yf) (4 - {x, yff (l - (x, yff {x, yf 

in X vanishes at each minimal vector but does not average to over the sphere of 
radius 2, because it is nonnegative on that sphere. 

The constant 2.50193 • 10~^ in Lemma [O can be made somewhat smaller for 
polynomials of total degree at most 8 (using the same proof), and that is the only 
case we need later. However, the present bound suffices. 



6.2. Intersection numbers. We will now use the fact that Ca is nearly a spherical 
10-design to determine the intersection numbers (still following the techniques of 
[DGS| ). As a computational aid, it is useful to know the following formula for 
averaging homogeneous polynomials over the sphere: 



1 



31(^23) 



523 



g{z)dz 



/R24 5(z)e I'^I'dz 



rdogse-'-' vol(523)r23 dr ' 



In the exact case of A24, we can compute the intersection numbers as follows. 
For each x,y £ C24 (the spherical code derived from the minimal vectors) with a 
specified inner product, we need to determine the number of z S C24 with specified 
inner products with x and y. Let P^{a,(3) denote this number when {x,y) = 7, 
(a;, z) = a, and (y, z) = (3. The cases 7 = ±1 simply amount to the valencies, which 
are determined automatically once the other intersection numbers are determined. 
For instance Pi(a,/3) = unless a = /3, and Pi (a, a) = ^^Po(a,/3) once we 
demonstrate that there every vector in C24 has a vector in C24 orthogonal to it 
(which follows from, say, showing that Pq,(0, 0) 7^ for each a 7^ ±1)- Hence we 
will focus on the remaining cases, i.e., 7 7^ ±1. 

For such a pair {x,y) with {x,y) = 7, a priori there are 49 unknowns P^{a,(3) 
fora,/3e {-1,-1/2,-1/4,0,1/4,1/2,1}. We first note thatP^(l,a) = P^ia^l) = 
Sa.-y, where 6 is the Kronecker delta. Similarly P-y(— l,a) = P-y(a, — 1) = 
Thus we can eliminate ±1 from consideration, which reduces the problem to finding 
only 25 unknowns. We will find 25 linear equations that determine these values. 
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Po(0,0) = 43164 
Fo(l/2,l/2) = 44 

p^/2(0,0) = 49896 
Pi/2(1/2,1/2) = 891 
Pi/2(l/2,-l/4) = 

Pi/4(0,0) = 44550 
Pi/4(l/2,l/2) = 275 
Pi/4(l/2,-l/4) = 275 



Po(0,l/2) = 2464 
Po(l/2, 1/4)- 1024 

Pi/2(0,1/2) = 891 
Pi/2(l/2,-l/2) = l 
Pi/2(l/4, 1/4) ^ 20736 

Pi/4 (0,1/2) = 2025 
Pi/4(l/2,-l/2) = 
Pi/4(l/4, 1/4) = 15400 



Po (0,1/4) = 22528 
Po(l/4, 1/4) = 11264 

A/2 (0,1/4) = 20736 
Pi/2(l/2, 1/4) = 2816 
^'1/2(1/4,-1/4) = 2816 

Pi/4(0, 1/4) = 22275 
Pi/4(l/2, 1/4) = 2025 
Pi/4(l/4, -1/4) = 7128 



Table 1. Intersection numbers for the Leech lattice minimal vectors. 



Consider the polynomials gi,j{z) = {z,xY{z,y)^ for i, j E {0, 1, . . . ,4}. Let S be 
the set {-1,-1/2,-1/4,0, 1/4, 1/2,1}. We then know that 

^ , ^, 196560 f , , , 

because C24 is a spherical 11-design (although even an 8-design would suffice). 
These equations can be solved to yield the unknown values of P-y{a,(3). Note that 
the right hand side does not depend on the choice of x and y, only on 7 = {x, y). 
Therefore we see that the solutions of these equations, which are the intersection 
numbers, are independent of x, y and only depend on 7. We solve one such system 
of equations for each value of 7, and the values of the intersection numbers are 
tabulated in Table [1] which is a complete list modulo the symmetries 

P^(a, P) = P^(/3, a) = P^(-a, -/3) = P-^{a, -/3). 

(These numbers are of course known, but they are not tabulated in standard refer- 
ences such as |CS2| . so we record them here for convenience.) 

What happens when we are not necessarily dealing with exactly A24, but rather 
with A? Suppose we have x,y € Ca and we want to determine the number of 
z e Ca with specified approximate inner products with them. Let Py(a,/3) denote 
the intersection numbers for Ca (which may depend on x and j/); here we use a, 
P, and 7 to denote exact elements of {0, ±1/4, ±1/2, ±1}, and the inner products 
from Ca are required to be approximately equal to them. Then we have 

weCA 
weCA 

196560 , , , , 

gij(z) dz 

23\ ^iji'y)^ 



VOl(523) 

196560 
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where 

G.,,(7)- / {z,ur{z,vydz 

with {u,v) = 7 (recall that {x,y) ~ 7). 

If these approximations are close enough, we should get about the same values 
for P-y(a, P) as we did for P-y{a, (3), and this will show that the intersection numbers 
must be the same. 

Lemma 6.3. Leta,(3e {0, ±1/4, ±1/2, ±1} and a, 5 e [-(1/2 + (t), l/2 + cr]U{±l} 
with max{|a — a\,\b — P\} < a < 0.1. Then for i, j > 0, 

\a'b^ -a'f3^ < (l + 2o-)cr. 

Proof. We begin by bounding — a*|. For i = or a G {±1}, a* — a% so we can 
assume that i > and a e {0, ±1/4, ±1/2}. Then 

i-l 

k=Q 

If we apply the triangle inequality (together with \a — a\<a and |q;|, \a\ < 1/2 + cr), 
we find that 

(6.1) |a^-a^| <i(1/2 + <t)'-V. 

The function n(l/2 + cr)"""'^ on nonnegative integers attains its maximum value 
when n = 2, from which it follows that |a* — a'| < (1 + 2a)a-. Note that this is the 
special case of the lemma in which P — b — ±1. Thus, we can henceforth assume 
that neither P nor 6 is ±1, and by symmetry we can assume the same for a and a. 
For the general case, we notice that 

a'P^ - a'V = {a' - a')P^ + a\P^ - V), 

so 

\a'P' - a'V\ < \a' - a'\ ■ + \a'\ ■ - V\. 
It now follows from (|6.ip that 

\a'p^ - a'V\ < (i + j)(l/2 + cr)'+^'-V. 
The right hand side is at most (1 + 2CT)a, as before. □ 
Now we need to make precise the errors in all three approximations in 
^ a'/3^"P^(a,/3)« ^ a.,,(w) 



196560 



5i J (z) dz 



The first follows from Lemma l6. 3 



a,f3eS ujGCa 



vol(523) y^,3 
196560 ^ 

V0l(^23)^^jiT) 



< ^ (1 + 2(t)o- = 196560(1 + 2o-)cr, 

wECa 



because \a'^P^ — a^b^\ < (1 + 2a)a for {w, x) = a w a and {w, y) —b~ p. 



16 



HENRY COHN AND ABHINAV KUMAR 



The second we have ah^eady estimated sufficiently weh in Lemma 16.21 in terms 
of |5ij|2- Because \gi,j{z)\ < 1 for aU z, we have 



It foUows that 



TT 



12 



12! V19958400 



^ 196560 f 

wECa ^ ' ^ 



< ■2. 50193-10'^ < 5.3841-10" 

V19958400 



Finally, one can check by straightforward computation of each case that if (x, y) 
differs from one of —1/2,-1/4,0,1/4,1/2 by at most a (where < ct < 1), and 
< «,j < 4, then 

196560 f , , , 

differs by at most 8190<t from what it would be if a were zero. (In fact, if one 
expands this quantity as a power series in cr, then the sum of the absolute values 
of the coefficients is at most 8190.) Therefore the error in the last approximation 
is at most 8190(t. 
Thus 

Ei r,'i ?, / r,^ 196560 „ , , „ 



where 



I A J I < 196560(1 + 2a)a + 8190cr + 5.3841 • 10-*^ < 6.7023 • 10"^ 



As before the values of P,y{±.\,a) and P^(a,±l) are known and are the same as 
the corresponding P^{±.l,a) and P-y(a,±l). Thus they serve as constants in the 
equation and do not contribute to the error. 

Let A be the matrix of coefficients for these equations. One can check that 
|A~^|oo = 7225. (Here, | • |oo denotes the oo-norm on matrices, which is induced by 
the too norm on vectors. It is the maximum over all rows of the sum of the absolute 
values of the elements in that row.) It follows that the intersection numbers in Ca 
differ by at most 

7225 • 6.7023 • 10"*^ < 0.05 

from those in C24. Because they must be integers, this proves that they are the 
same as in C24 (in particular, they do not depend on the choice of x and y, although 
what we have proved is far stronger). 

The computer file verifyrest.txt carries out all these calculations (as well as 
those from several other points in this paper). In it, we assume without loss of 
generality that 7 > because P^{a,(3) = P^^{a, —f3). 

6.3. Uniqueness of association scheme. The spherical code C24 determines a 6- 
class association scheme A24: if we partition elements {x, y) of C24 x C24 with x ^ y 
according to their inner products {x,y). We can similarly form the association 
scheme of the spherical code Ca coming from A, where this time we group 
elements according to their approximate inner products. We wish to show that 
these two association schemes are isomorphic. We will need to know that this 
6-class association scheme ^24 is uniquely determined by its size, valencies, and 
intersection numbers. That can be proved as follows. 
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Let N = 196560, let C = 6/iV = 1/32760, and let ui,...,un be the minimal 
vectors of the Leech lattice. The following lemma restates the (known) fact that 
A24 is strongly eutactic. We provide a proof here to make this article more self- 
contained. 

Lemma 6.4. For every x G M?^ , 

N 

{x,x) = C y^^{x,Ui)'^. 

1=1 

Proof. Because the minimal vectors of A24 form a spherical 2-design, the polynomial 
y {x,y)'^ has the same average over {ui, . . . ,un} and the sphere of radius 2 in 
M^^, which is 4 times the average over the unit sphere. To average over the unit 
sphere, it will be convenient to work with orthonormal bases of R^*. For each 
orthonormal basis ei,...,e24, we have |xp = ^^(x, e^)^. If we average over all 
orthonormal bases, then each of ei, . . . , 624 is uniformly distributed over the unit 
sphere, and therefore the average oi y 1-^ {x,y)'^ over the unit sphere is \x\'^/2A. It 
follows that the average over the sphere of radius 2 is |a;p/6, so 

1 ^ 

-Y.{a:,u,y = \x\y6, 
1=1 

as desired. □ 

Theorem 6.5. There is only one 6-class association scheme with the same size, 
valencies and intersection numbers as the association scheme of minimal vectors of 
the Leech lattice. 

Proof. Let A = {a^ : 1 < i < 196560} be such an association scheme, with 
partitioned into classes as 

^AiD A/2 u Ai/i u A u Ai/4 u Ai/2 u A-i 

(labeled according to the corresponding inner products in the Leech case, when the 
minimal vectors are rescaled to lie on the unit sphere). 

Let Ai, Ai/2, ■ • ■ , A-i denote the adjacency matrices of the identity relation and 
the six classes; in other words. 



(^0)4,; 



1 if {ai,aj) e Aa, and 
otherwise. 



These matrices are symmetric and commute with each other. Their span forms 
an algebra, called the Bose-Mesner algebra, whose product is determined by the 
valencies and intersection numbers. (See [BH} p. 755].) Namely, 

AaAfj = ^P^(a,^)A■ 
7 

Furthermore, note that because the Aq,'s have only or 1 as entries, are not identi- 
cally zero, and sum to the all I's matrix, they must be linearly independent. Thus, 
the Bose-Mesner algebra is completely determined by the valencies and intersection 
numbers, with no additional relations possible. 
Let 

P = 4C Ml + iAi/2 + JA/4 - J^-1/4 - ^^-1/2 - ^-1 
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If the association scheme A comes from the Leech lattice, then P is C times the 
Gram matrix of the 196560 minimal vectors (not rescaled to the unit sphere). We 
claim that P is a projection matrix, i.e., P^ = P. One way to check that is to use 
the valencies and intersection numbers to compute P^ and verify that it equals P. 
That is somewhat cumbersome to check by hand, so we will instead give a longer 
but conceptually simpler proof. 

First, we check it in the case of the actual Leech lattice association scheme, using 
Lemma 16.41 From 

N 

it follows that 

{x, y) ^]^{{x + y,x + y) ~ (a;, x) - (y, y)) 
C ^ 

i=l 
N 

= C^{x,u^){y,Ui). 
1=1 

Therefore the entry of P^ is 

N 

(P^).,.j = C^^(Ui,Mfc)(Wj,Mfc) C{u.i,Uj) = Pij, 
k=l 

SO for the Leech lattice P is a projection matrix. Because the structure of the Bose- 
Mesner algebra is completely determined by the valencies and intersection numbers, 
the same must always be true. 

Thus, P is always a projection matrix (in fact, an orthogonal projection because 
P is symmetric) . The trace of P is AC times the trace of Ai , because no other A^ 's 
have entries on the diagonal. The trace is therefore 4iVC = 24, so P projects onto 
a 24-dimensional subspace. 

Consider the images of the N = 196560 unit vectors (namely ei, . . . , eigeseo) 
under P. Their inner products are simply the entries of P, since {Pei,Pej) = 
{ei,P^ej) = {ei,Pej) = Pij, and if one rescales the vectors by l/ViC so that 
they lie on the unit sphere, then the result is a 196560-point kissing configuration 
in R^**. (Appendix IB] reviews the kissing problem.) The only such configuration 
is the kissing configuration of the Leech lattice (see [BSj . which was reprinted as 
Chapter 14 of [CS2j ). up to orthogonal transformations of M^*. Thus, P/C must 
be the Gram matrix of the minimal vectors in the Leech lattice. 

It follows that A is isomorphic to A24 (in particular, ai t-^ P{ei)/ VaC yields 
an isomorphism). Thus, ,424 is determined by its size and intersection numbers, as 
desired. □ 

7. Inner product bounds 

We will first use the intersection numbers and isomorphism of association schemes 
to prove better bounds on a. We already know that for all nearly minimal vectors 
u, 

4 < {u,u) < 4(l + £)2 < 4 + 9e. 
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Lemma 7.1. For nearly minimal vectors u,v with {u,v) ~ 2, we have 

2 - 5£ < (m, u) < 2 + 9e. 
Proof. We know that u — v is nearly minimal, so 

{u- v,U'- v) < 4(1 + e)^. 

It follows that 

2{u, v) > -4(1 + 2e + e^) + {u, u) + (i;, v) 

> 8-4(l + 2£ + e2) 

> 4 - lOe, 

which gives us one of the inequalities. Similarly, 

{u — v,u — v) > 4 

gives us 

2{u, w) < -4 + (u, u) + {v, v) 

< -4 + 8(l + 2£ + £2) 

< 4+ 18£, 

which is the other inequality. □ 

Note that these inequalities could be made slightly sharper, but we prefer simpler 
numbers. 

Lemma 7.2. For nearly minimal vectors u, v with {u, v) w 0, we have 

-lAe < (u, v) < 14£. 

Proof. Since Po(l/2, 1/2) = 44 ^ 0, we can find w with {u, w) fa 2 and (v, w) « 2. 
Then v — w is nearly minimal and we know from the previous lemma that 

2 - 5£ < (u, w) < 2 + 9£ 

and 

-2 - 9£ < {u, v-w) < -2 + 5£. 
Adding these inequalities gives us the result. □ 

Lemma 7.3. If u,v e A24 are minimal vectors with {u,v) = 1, then there are 
minimal vectors Wi, W2, andw^ satisfying {u,Wi) — 2, {v,Wi) — 0, and {wi,Wj) = 
fori ^ j. 

Proof. We will use the fact that the automorphism group Coq of the Leech lattice 
acts transitively on pairs of minimal vectors with a fixed inner product between 
them (see Theorem 3.13 in [T] for a proof). Thus it suffices to consider the case of 
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1 


u — 




V 


1 






1 


Wl = 


71 




1 


W2 


71 




1 


W3 = 


71 



a particular pair of minimal vectors with inner product 1. Let 

= (1,1,. ..,1,1,-3), 

= (0,0,..., 0,-4, -4), 

= (2,2,2,2,2,2,2,2,0,0,. ..,0), 

= (0,0, ...,0,4,-4), and 

I 

= (0, 0, 0, 0, 0, 0, 0, 0, 2, 2, 0, 0, 2, 2, 0, 0, 2, 2, 0, 0, 2, 2, 0, 0). 

(see |CS21 p. 131] for a description of the minimal vectors). It is easily checked that 
the inner products are as desired. □ 

Lemma 7.4. For nearly minimal vectors u,v with {u,v) ~ 1, we have 

1 - 72e < (u, v) <1 + 75e. 

Proof. We know that the association schemes of the Leech lattice A24 and our 
given lattice A are the same. Let u',v' be the corresponding vectors in the Leech 
lattice (corresponding via some fixed isomorphism of association schemes). Since 
(u', v') = 1, we know that we can find w'l, W2, and ^3 in the Leech lattice satisfying 
{u',w[) = 2, {v' ,w[) = 0, and {w[,w'^) = by Lemma Let wi , ^2 , ^3 be 

the corresponding vectors in A. Then the relations (u, wt) ~ 2, {v, Wi) ~ 0, and 
(wi, Wj) ~ must hold in A. It follows by a short computation that 2u — u — wi — 
W2 — W3 is a nearly minimal vector. Therefore 

4 < (2u — V — Wl — W2 — Wz,2u — V — Wl — W2 ~ W3) 

= 4(u, u) + {v, v) + ^(wi, Wi) - 4(u, v) 

i 

- 4 ^(m, w,) + 2 ^(v, w,) + 2 ^(w,, u.j). 

i i i<j 

It follows that 

4(u, w) < -4 + 4(u, u) + {v, v) + ^(wi, Wi) 

i 

-4^(u,w,) +2^(t>,w,) +2^(w,,Wj) 

i i i<j 

< -4 + 8(4 + 9e) - 12(2 - 5e) + 12(14e). 
Thus, (u, v) < I + 75e. Similarly, 

4 + 9£ > (2m — V — Wl — W2 — W3,2u — V — Wl — W2 — w^) 
= 4{u, u) + (w, v) + ^(wj, Wi) - 4(m, v) 



4^(u,Wi) +2^(^;,Wi) + 2^(wi,Wj) 



i<j 
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It follows that 

4{u, v) > -4-9e + 4(it, u) + {v, v) + ^(wi, Wi) 

i 

i i i<j 

> -4 - 9£ + 8(4) - 12(2 + 9e) + 12(-14£). 

Thus, {u, v)>l- (285/4)e > 1 - 72e. □ 

Wc have proved the following proposition. 

Proposition 7.5. If u,v € A are nearly minimal vectors, then {u, v) differs from 
an element 0/ {0, ±1, ±2, ±4} by at most 75e. 

8. A BASIS OF NEARLY MINIMAL VECTORS 

We wish to prove that A must have a basis of nearly minimal vectors. We first 
prove that the nearly minimal vectors span A. Let x S A be as small a vector 
as possible without being in the span of the nearly minimal vectors. Then there 
does not exist a nearly minimal vector u such that |u — a;| < |a;|. We know that 
|a;| > ^6(1 — m) ^iid 2 < |m| < 2(1 + e) for each nearly minimal u. Because 
I'M — a;| > \x\, we have 

{%x) < |m|V2 < 2{l+sf. 

Consider the unit vectors x/\x\ and u/\u\. We have 

u X \ 2(1 + e)2 1 

H' R/ " %/6(1-m)-2 ^ 2' 

If we extend the spherical code Ca of vectors of the form u/\u\ to Ca U {a;/|a;|}, then 
it will contain 196561 vectors without changing the minimal angle, and we have 
seen that that is impossible. Thus, the nearly minimal vectors do span the lattice. 

The same argument (with e = n = 0) also proves the following fact: if the kissing 
configuration of a lattice is an optimal kissing configuration for its dimension, then 
the lattice is spanned by its minimal vectors. 

Now let B be l/VS times the matrix 

/ 4 -4 0000000000000000000000\ 
' 440000000000000000000000' 
4 4 
400400000000000000000000 
400040000000000000000000 
400004000000000000000000 
400000400000000000000000 
222222220 000 0000 
400000004000000000000000 
400000000400000000000000 
400000000040000000000000 
222200002222000000000000 
400000000000400000000000 
220022002200220000000000 
202020202020202000000000 
20 2200220022 02 00000 
400000000000000040000000 
202020022200000022000000 
200222002020000020200000 
22 20202002000020020000 
2 2 2 2 2 2 2 2 
000000002200220022002200 
000000002020202020202020 
-3 11111111111111111111111 
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The rows of B form a basis of A24 consisting of minimal vectors (see Figure 4.12 
in [CS2j ). One can compute the inverse matrix and check that its entries are 
integers divided by \/8- The largest entry of B^^, in absolute value, is — 13/\/8, so 
|B-i|oo < 24 -13/78. 

Let u'l, . . . , U24 be the rows of B. Every minimal vector u' is a linear combination 
Ciu'j^ of this basis, and the coefficients are bounded by 



Consider the corresponding nearly minimal vectors ui, . . . , tt24 of A (under an iso- 
morphism of association schemes). We next prove that they form a basis. 

We need only check that ui, . . . , U24 span the nearly minimal vectors of A. If 
u is any nearly minimal vector, the isomorphism of association schemes gives us 
numbers Ci, . . . , C24 such that u should be J^i '^i'^^i (^-^-i the corresponding equality 
is true for the Leech lattice). We first check that c^u^ is a nearly minimal vector, 
and then that it equals u. To check that it is nearly minimal, we just need to know 
that all inner products of nearly minimal vectors are within 75£ of what they are in 
the Leech lattice (Proposition 17. 5p . When we compute the norm of CiUi, each 
inner product (u^, Uj) could be off by as much as 75e, and multiplied by up to 156^. 
There are 24^ such pairs, for a total error of at most 156^ • 24^ • 75e = 1051315200e, 
which is minuscule (less than 10~^^). Because the next smallest vectors beyond the 
nearly minimal vectors have norms at least •\/6(l — m); conclude that '^^ciUi 
must be nearly minimal. 

To check that u = CiU,;, we need only compute their inner product and verify 
that it is approximately 4. This time, the maximum error is 24 • 156 • 75£, which is 
again small enough by a huge margin. 

Thus, Ml, ... , M24 form a basis of A, and their inner products are all within 75e 
of what they would be in the Leech lattice. 



In this section, we will prove in detail that the Leech lattice is locally optimal, 
and provide quantitative bounds. We will follow the notation and techniques from 
|GLj closely. The basic result is Voronoi's theorem, which says that a lattice is a 
local maximum for sphere packing density if and only if it is perfect and eutactic. 
In this section of the paper only, A will denote an arbitrary lattice in R" . 

It is convenient to work in terms of the quadratic form Q associated to A. Choose 
a lattice basis {h}, and for a vector a; G K" (with coordinates xi, . . . , Xn) define 



The matrix of Q is 5 = {si,jYi^j^i, where Sij = {bi,bj). 

Let M denote the minimal nonzero norm of A, and let ui, . . . ,un G be the 
coefficient vectors of the minimal vectors in terms of the basis {bi}. Thus, for 
l<i<N, 



Recall that Q is perfect if these equations completely determine the quadratic form 
Q. Equivalently, every quadratic form that vanishes at ui,...,ujv must vanish 
everywhere. 



|c.| < \B-^\oo\u'\oo < (24 • 13/VS) ■ (4/V8) = 156. 



9. Local optimality of A. 



-24 
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Let D denote the determinant of S, and let S — denote the adjoint 

matrix, where 

, _ dD 
''^ ~ 9s,j ' 

Strictly speaking this is an abuse of notation, but of course it means we take the 
partial derivatives of D as if the entries of S were variables, and then substitute 
their actual values. 

In other words, SS is D times the identity matrix. (It might seem that a trans- 
pose is missing, but note that all our matrices are symmetric.) Let Q be the 
quadratic form with matrix 5'. Then A is eutactic if there are positive numbers 
di, . . . ,dN such that for all a; S R" , 

N 

Q{x) = ^ dk{uk,x)^. 
fc=i 

It is known that A24 is perfect and eutactic, with di = ■ ■ ■ = digeseo = 1/32760. 
For completeness we sketch a proof. Lemma 16.41 proves that the Leech lattice is 
eutactic. 

Lemma 9.1. The Leech lattice A24 is perfect. 

Proof. Suppose Q is a quadratic form that vanishes on the minimal vectors. Then 
the symmetric bilinear form B corresponding to Q satisfies B(ui,Ui) — 0. We first 
show that B{u, v) = for all minimal vectors u, v. If {u, v) — 2 we use the fact that 
u — V is minimal to see that 

B{u, = ^ {B{u, u) + B{v, v) - B{u - v,u~ v)) ^ 0. 

If {u,v) — 0, then since Po(l/2,— 1/2) = 44 7^ 0, we can find w with {u,w) — 2 and 
{v,w) — —2. Then t; + w is nearly minimal and we know 

— B{u, w) — Biu, V + w)] 

subtracting these two gives us the result. If = 1, then by Lemma 17.31 there 

are minimal vectors wi, W2, and 1^3 with {u,Wi) = 2, {v,Wi) = 0, and {wi,Wj) = 
for i j. It follows that 2u — v — wi — W2 — W3 is a minimal vector. Therefore 

= B{2u — V — wi — W2 ~ W3, 2u — V — wi — W2 — w^) 
= 4B(u, u) + B{v, w) + ^ B{wi,Wi) - 4:B{u, v) 

i 

i i i<j 

= -4:B{u,v)+0. 

This forces B{u, v) — 0. Finally, the result for (u, v) < follows from the above by 
noticing B{u,v) = —B{u, —v) and (m, —v) > 0. 

To conclude the proof, we use the above information on a basis of minimal vectors 
to see that B, and hence Q, is identically zero. □ 

We begin by proving one direction of Voronoi's theorem. This proof is the one 
given in [GLi §39], where one can also find a proof of the converse. We give the 
details of this direction because we will need to examine it in detail to derive 
quantitative estimates, and because it is the only direction we will need. 
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Theorem 9.2 (Voronoi | Vo| ) . //A is perfect and eutactic, then it is a strict local 
maximum for density. 

Proof. We wish to show that if A is perturbed shghtly (other than simply by scahng 
and isometries), then D~^/"Af must strictly decrease. We perturb S = (sij) by 
changing it to (s^j + ptij), where tj^i = tij and p > is small. We assume (tij) 
is not identically zero. Let Qp denote the corresponding quadratic form, and let 
Dp be the determinant of the matrix {sij + ptij). We will show that for any fixed 
matrix (i^.j) not proportional to S, if p is sufficiently small, then there exists a k 
such that 

Because Qp{uk) is no smaller than the minimal norm of Qp (i.e., the minimum 
value of Qp on Z" \ {0}), this inequality is what we want. 

First, note that without loss of generality we can assume that 

(9.1) J2i,^jt,^^^0. 
The reason is that 

Every perturbation can be broken up into the sum of a perturbation proportional to 
S and a perturbation satisfying ()9.1|) . If we deal with the latter, then we can ignore 
the former (which rescales the lattice but does not change its packing density). 
Thus, we assume ()9.ip from now on. Then the determinant Dp is given by 

Dp = det(sjj + pti^j) = det(sij) + p^iijSij + 0[p^) 

= D + 0{p^) as p ^ 0. 
Because A is eutactic, 

JV 

k=l 

for all X e ffi.". The associated symmetric bilinear form is 

N 

^dk{uk,x){uk,y), 
fc=i 

from which it follows that 

N 

Si,j = ^ dk{uk)i{uk)j, 

k=l 

where {uk)i denotes the i-th coefficient of Uk. Therefore 

N 

dk^^tij{uk)i{uk)j = ^ Sijti^j = 0. 

k—l ij i,j 

Because A is perfect, the inner sum on the left hand side in the above equation 
cannot vanish for all k. Therefore, there exists a k for which it is negative, say 

^tij{uk)i{uk)j < -a 
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with a > 0. Then 

Qp{uk) = '^{si,] + pU,j){uk)i{uk)] < M - pa, 

and hence 

i/"QpK) < i:»-i/"Af (1 - pa/M)(l + 0(p2))-i/« < D-^I'^M 

if p is positive and smaU enough. This proves that A is a strict local optimum for 
density when (sij) is perturbed in the direction of {tij). In fact, the choices of 
a and the implicit constant in the big-O can be made uniformly in (tij), given 
Si j l^jj P — 1- Thus A is a strict local optimum for density. □ 

We now compute numerical bounds on the perturbations. We use the same basis 
of A24 as before. The corresponding Gram matrix is 



v 



Suppose we assume that the Gram matrix entries are perturbed by pU j where 
max{|ii,j|} < 1 (and ()9.ip holds — we will eventually deal with the case where it 
does not). Then, using the same notation as in the proof of Theorem [ 

Dp-D = J2 p'^™^ det(T) det(f ), 



where T ranges over all nonvacuous minors of (tij) and det(T) is the corresponding 
cofactor of the matrix S. (The D term corresponds to the case where T is vacuous, 
i.e., contains no rows or columns.) This expansion follows from combining the 
Laplace expansion (see jMu|, §95]) with multilinearity, and is known as Albeggiani's 
theorem [Mul §96]. 

Now let dimT = k. Then |det(T)| < fc''/^ by Hadamard's inequahty, since the 
length of each row of T is at most ^/k (see (7.8.2) in |H J| ) . It follows that the 
absolute value of the sum of p'' det(T) det(f') over fc-dimensional T is bounded by 
k^/'^p^Ak, where Ak is the sum of the absolute values of the (24 — fc)-dimensional 
minors of the Gram matrix. For fc > 3 we use the simple bound Ak < i^^) ■ (4^ + 
(24-fc-l)(22))(24-fc)/2 (the first factor is the number of (24 -A:) x (24 -A:) minors, 
and the second is the above bound on the determinant of the cofactor, using the 
fact that the largest entry in each row of 5 is 4 and the other entries are at most 
2 in absolute value). For k — 2 we explicitly compute the sum of absolute values 
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of the 22 X 22 minors, and find that it is 818153. Putting aU this together, we see 
that for < p < lO^^o^ 

(9.2) Dp> l-p2(22/2.818153 + 2-10*), 

where the 2 • 10® term bounds the contribution from all higher powers of p. (Recall 
that D ~ 1 for the Leech lattice.) See the computer file verifygram.txt for the 
details of this calculation. 

Next, we find an a that works for any choice of tij such that max{|fij|} = 1. 
This is done by linear programming as follows. 

We find an a > such that for all («0i Jo); and for all tij subject to the constraints 
-1 < tuj < 1 for ^ {io,jo), ti,] = tj.i, tio,]o = 1' and Y.i,j = 0' the 

following inequality holds for some minimal vector Uk'. 



E 



We do the same for tip, jo = — 1. All these linear programs can be solved by com- 
puter, and it appears that a — 1/23 works and is the largest possible value of a. 
However, that is the result of floating point calculations that are not rigorous. We 
will be content with proving (without computer assistance) that a — 4/1055 satis- 
fies these properties. This weaker bound will suffice for our purposes and is proved 
in Section [TUl 

We conclude that for every perturbation by tij where max{|tij|} = p and 



J2i SijUj = 0, there exists a k such that for < p < 10 



(9-3) ^C'/?/.^M^ (1 - '^/1055)(1 - (2 . 10« + 2 . 818153)pY'^'\ 



We used here that M = 4 for A24. The upper bound in (19. 3p is strictly less than 1 
when < p < 10^^". (Note that for notational convenience we have absorbed the 
factor p into the perturbations iij .) 

The last remaining issue is that our perturbation may not satisfy j — 0- 

Suppose our perturbed matrix entries are Si^ + r]ij . Let 

A — ^ijVij ■ 

By Proposition 17.51 \rji^j\ < 75e. It follows that |A| < 152100e, since the sum of 
the absolute values of the entries of S is 2028. 

If we divide the quadratic form by l+A/24, which is nonzero, then it is equivalent 
to A24 perturbed by 

Si.j + rjij 



l-f A/24 



where now 



because 



^s.js^^j = 24D = 24. 



To conclude our proof, we need only check that \tij\ < 10~^°. (Note that if tij = 
for all then the perturbed quadratic form is proportional to the original one 
and therefore equal to it because they have the same determinant.) 
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Using |A| < 152100e, \r]ij \ < 75e, and |sij| < 4, we have 
77,; J - s»,iA/24 



(9.4) 



75£ + 4- 152100e/24 _ _ ,3 

< — < 1.8 • 10"^^. 

- 1 - 152100e/24 



1 + A/24 

Because 1.8 • 10^^^ < 10^^", wc find that the dense lattice A is close enough to A24 
to conclude that A is either the same as A24 (up to isometrics of R^"*) or strictly 
less dense than A24. 

This completes the proof of our main theorem, except for the postponed com- 
putation of a in Section [TOl 

Theorem 9.3. The Leech lattice is the unique densest lattice in M^^, up to scaling 
and isometrics o/R^^. 

Note that the reason why scaling ambiguity appears in the theorem statement 
but not in the above proof is that we fixed |A| = 1. 

10. Computation of a 

Suppose zq, Jo G {1, 2, ... , 24} and t = ±1. We wish to find a number a > 
such that whenever —1 < tij < 1 for {i,j) ^ (*0j jo), ti.j = tj,i, tio,jo = ^-i^d 
Si j ^i,3^i,3 ~ 0, there is a fc such that 



'^U,]{uk)i{uk)j < -a. 



We will see that we can take a = 4/1055, whatever zo, jo, and t are. 

Lemma 10.1. Each minimal vector of the Leech lattice is contained in a set of 
twenty-four orthogonal minimal vectors. 

In fact, more is true: the minimal vectors can be partitioned into 4095 sets 
{±vi, . . . , ±i;24} with Vi and Vj orthogonal for i ^ j. See footnote 3 of |E1[ p. 6] 
for an elegant proof. 

Proof. Because the automorphism group of the Leech lattice acts transitively on 
the minimal vectors, we need only verify that there exists a set of twenty-four 
orthogonal minimal vectors. Let vui be the vector 

-^(0,. .. ,0,4,4, 0,...,0), 

where only coordinates i and i + 1 are nonzero, and let Vi be the vector 

-^(0,. ..,0,4, -4,0,. ..,0). 

These are all minimal vectors in the Leech lattice, and wi, vi, W3, . . . , ^23, V23 is 
an orthogonal basis of R^"' . □ 

Let T denote the matrix {tij). We will also write T{v) = {v,Tv) and T{u,v) = 
{u, Tv). In these terms, our goal is to show that given the assumptions on T, there 
is some k such that T{uk) < — a. 

We will make use of the following lemma, which depends only on the hypotheses 
listed in its statement: 
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Lemma 10.2. Let wi, . . . , W24 be orthogonal minimal vectors in the Leech lattice. 
VJ2t,j Si,jtij = 0, then 



24 

E 

2=1 



T{v,) = 0. 



Proof. Let B be the matrix whose rows are the coordinates of ui , . . . , U24 relative 
to the basis we have chosen for the Leech lattice. Then we see by orthogonality 
that BSB* = 4/. Thus, S = AB'^iB*)-^. It follows that 



Tt{BTB^) = Tr(S*BT) = 4Tr(5"^T) = 4 S^jt^^ = 0, 
which implies 



24 

1=1 



□ 



Let us rephrase our basic problem, and slightly weaken the hypotheses, as follows. 
Suppose J Sijtij = 0, and we are given minimal vectors x, y such that 

=/3e {±4,±2,±1,±0} 

and T{x,y) = t ^ 0. We wish to find a > such that there must always exist a 
minimal vector w with T{w,'w) < —ct\t\ under these hypotheses. We will apply it 
with t = ±1. (Note that we are no longer assuming \ti,j \ < 1.) 

If we prove a bound for a certain t) we automatically get the same bound of 
a for the case of {—(3, —t) (just consider the pair [x, —y) instead of (a;, y)). Thus it 
suffices to prove bounds of a for the cases when /3 is nonnegative. 

Lemma 10.3. If f3 — 4 and t <0 then we can take a — 1. 

Proof. The hypothesis says that x — y and T{x,x) = t < so we just take w = 

X. □ 

Lemma 10.4. If (3 ~ A and t > 0, then we can take a = 1/23. 

Proof. The hypothesis says that x = y and T{x,x) — t > 0. By Lemma llO.ll 
the Leech lattice contains orthogonal minimal vectors Vi, . . . ,V24 with vi — x. 
Then Lemma [10.21 together with T{vi) = t, implies that T{vi) < — 1/23 for some 
{2,...,24}. □ 

Lemma 10.5. If (3 = 2 and t <Q, then we can take a = 2/25. 

Proof. Let x^y be minimal vectors such that {x^y) — 2 and T{x,y) — t < 0. We 
know that a; — y is a minimal vector. Also 

T{x - y) - T{x, x) + T{y, y) - 2T{x, y) = T{x, x) + T{y, y) - 2t. 

Thus either T{x,x) < 2t/2b = -2|t|/25 or T{y,y) < 2t/25 or T{x -y,x-y) > 
At/25 - 2t = -46i/25 = 46|i|/25. However, in the last case we see by Lemma fTO^ 
that T{v,v) < -(l/23)(46|t|/25) = -2|t|/25 for some minimal vector v. □ 

Lemma 10.6. If (3 — 2 and t > 0, then we can take a = 2/47. 
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Proof. Let x, y be minimal vectors such that (x, y) =2 and T[x, y) = t > 0. Again 
we have 

T{x -y,x-y) = T{x, x) + T{y, y) - 2t. 

Thus either T{x — y,x — y) < —2t/A7 or one of T{x,x) and T{y,y) is at least 
{2t — 2t/47)/2 = 46^/47. Then by Lemma [10.41 we see that there exists a minimal 
vector V with T{v,v) < {~l/23){m/47) = -2\t\/47. □ 

Lemma 10.7. If /3 = Q, then we can take a — 1/60. 

Proof. By possibly exchanging {(3,t) = (0,i) with (0, — i) we can assume t > 0. 
Let x,y be minimal vectors such that {x,y) = and T{x,y) = i > 0. We have 
computed the intersection numbers for the Leech lattice. The intersection number 
Po(l/2, 1/2) is 44, so there exists a minimal vector w with {x, w) ~ 2 and {y, w) = 2. 
We compute that 

{x + y - w,x + y - w) = {x,x) + {y, y) + (w, w) + 2{x, y) - 2{x, w) - 2{y, w) 
=4+4+4+0-4-4 
= 4, 

so a; + 2/ — w is a minimal vector. Then we compute 

T{x + y — w) — T{x, x) + T{y, y) + T{w, w) + 2T[x, y) — 2T[x, w) — 2T{y, w) 

= T{x, x) + T{y, y) + T{w, w) + 2t - 2T{x, w) - 2T{y, w). 

If T(a;, x) or T{y, y) or T{w, w) is at most —t/60 we are done. Similarly if T(x + y — 
w, x+y — w) > 23i/60 we are done, by Lemma ll0.4[ so we may assume none of these 
is true. Then we see that T{x,w) + T{y,w) > 1/2(2^-3^/60 - 23^/60) = 47^/60. 
It follows that one of the summands is at least 47t/120, say T{x,'w) without loss 
of generality. But since {x,w) = 2, an application of Lemma 110.61 finishes the 
proof. □ 

Lemma 10.8. If (3 = 1 and t > 0, then we can take a — 4/1055. 

Proof. Let x and y be minimal vectors with {x,y) = 1. By Lemma 17.31 there are 
minimal vectors wi, W2, and W3 satisfying {x, Wi) = 0, {y, Wi) = 2, and {wi,Wj) = 
for i 7^ j. It follows that z — 2y ~ x ~ wi W2 ~ has norm 4, so it is a minimal 
vector. We have 

T{z, z) = 4T{y, y) + T{x, x) + ^ T{wi, Wi) - 4r(a;, y) 

■i 

x,Wi) + 2^T{w^,Wj). 

i i i<j 

Since we know T(x, y) — t, we get 

At = 4r(y, y) + T{x, x) +Y, T{w^, w^) - T{z, z) 

i 

- 4 ^ T{y, Wi) + 2Y, Tix, w,) + 2j2 T{w,, w^). 

i i i<j 

Now we will assume that T(v, v) > —at for all minimal vectors v, and prove that a 
cannot be less than 4/1055. It follows from this assumption and Lemma [10.41 that 
T{y,y) < 2'iat. Similarly T(x,x) and T{wi,Wi) are at most 23at, and —T{z,z) < 
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at by hypothesis. The inner products (y, Wi) are 2 so by Lemma llO.Sl —T{y, Wi) < 
25at/2. Finally, T{x, Wi) and T{wi, wj) are at most 60at by Lemma riO.TI Therefore 

4t < 8 • (23a<) + at + 3-4 - {25at/2) + 2 • 6 • {60at) = 1055ai, 

and hence a > 4/1055. It follows that T{v,v) < — 4i/1055 for some minimal vector 
V, since otherwise a could be decreased. □ 

Lemma 10.9. If (3 — 1 and t < 0, then we can take a — 4/1033. 
Proof. With the same notation as in the proof of Lemma FlO.SI we have 
T{z, z) = AT{y, y) + T{x, x)+J2 T{w^,w^) - 4T(a;, y) 

i 

~ AY^T{y,Wi) + 2Y^T{x,Wi) + 2Y,T{w,,w,). 

i i i<j 

Since we know T{x, y) — t < 0, we get as before 

-4t = -4:T{y, y) - T{x, x) T{wi,Wi) + T{z, z) 

i 

+ AY,T{y, w^) - 2 ^ T{x, w.) - 2 ^ T{w^, Wj). 

i i i<j 

We now use the same strategy as in Lemma [10.81 we assume that for all minimal 
vectors v, T{v,v) > ~a\t\ ~ at. It follows that —T{y,y) < a\t\ = —at, and the 
same holds for T{x, x) and T{wi, Wi). By Lemma [10.41 T{z, z) < —23at. This time 
by Lemma [10. 6[ T{y,Wi) is at most —i7at/2, whereas by Lemma [10. 7[ T{x,Wi) 
and T{wi,Wj) are at least QQat, so their negatives are at most — 60ai. Therefore 

< -{23at) - Sat + 3 -A - {-47at/2) + 2 • 6 • (-60at) = -1033ai, 

and hence a > 4/1033 after cancelling —t which is positive. □ 

We conclude that a = 4/1055 satisfies the properties we stated at the beginning 
of the section. 



11. The case of Es 

A very similar proof shows that the Eg lattice is the unique densest lattice 
packing in R^. Since the details of the proofs are analogous, and the result was 
already known, we merely sketch them in this section. 

In Eg, there are 240 minimal vectors of length \/2, if we normalize the length as 
usual so that the lattice is unimodular. Let A be a lattice of covolume 1 that is at 
least as dense as Eg,. As in the case of the Leech lattice, we find a suitable radial 
function / with r < \/2(l + 1.2 • 10"^^), which proves that no sphere packing in R 
can exceed the density of Eg by a factor of more than 1 + 10~^^. 

This function allows us to show that for 7^ |a;| < the length of x is restricted 
to the set 



[V2, V2(l + s)) U (2(1 - fi), 2(1 + fi)) U (V6(l - ly), V6(l + ly)), 
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where we have 

e = 1.45 • 10~l^ 

1.03- 10"^ and 
V = AAA - IQ-^. 

The details of this calculation are in the accompanying PARI file ESverifyf .txt. 
All the remaining calculations from this point on for the £"§ case are verified in the 
Maple file E8rest.txt. 

Define a nearly minimal vector to be a vector with length in [\/2, V2(l + e)]- 
We form a spherical code by rescaling the nearly minimal vectors to lie on the unit 
sphere S'^ . 

We use the polynomial 

= + (- + 0'-^ (--(i-2(TT7? 

to show as before that there are at most 240 nearly minimal vectors. Similarly, the 
analogue of Lemma [4.41 goes through with a different function defined on R^, and 
proves that there are exactly 240 nearly minimal vectors. 

11.1. Spherical code. The analogue of Section [5] is that all the inner products 
between the normalized nearly minimal vectors u/\u\ must be either ±1 or at most 
6 • 10~^ from some element of the set {—1/2,0, 1/2}. Then further analysis gives 
us the following better bounds. 
For {u, w) 1 we have 

2- (l + e)2 / u V \ 2 1 



2(1 + .)^ -\RH/-^'^^^ 2' 
whereas for (u, w) « we have 

1- (1 + Ai)2 / u w \ / , ,2 



We conclude that a < 8.89-10-^, where (7 is the maximal error in the inner products 
from the spherical code. 



11.2. Intersection numbers. The analogue of Lemma [6.21 is easily shown using 
the polynomial fg. 

Lemma 11.1. If g : S"^ M. is a polynomial of total degree at most 6, then 



Y] 9{z) / 9{z)dz 

where \g\2 denotes the norm on l?(S'^"). 



< 3.48- lO-^'Igla 



Now, since the kissing configuration of Eg, is a 7-spherical design, we find that 
the intersection numbers of Eg, can be obtained by solving the linear system of 
equations 

J2 a'P'P^{a,P) = / 9^,J{z)dz. 

a,/3eS ^ ' 
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for i,j g {0,1,2}, where S = {0,±1/2,±1} is the set of possible inner products. 
As before, we know the values of P^{zLl,a) and P^{a,ztl). When we perform the 
same calculation for A, the value of 



differs by at most 30(t from the corresponding value for Eg. If we apply Lemma l6.3l 
and compute the error introduced into the system of equations by going from 
to A, we get a bound of 



The oo-norm of the inverse matrix is 100, from which it follows that the intersection 
numbers in A differ from those in Eg by at most 0.44. Because that number is less 
than 1, they must be the same. 

11.3. Association scheme. The proof of uniqueness for the Leech lattice associ- 
ation scheme depended only on the eutaxy of the Leech lattice and the uniqueness 
of its kissing arrangement. The same hold for the Eg lattice. With N = 240 and 
C — A/N — 1/60 we have the following lemma. 

Lemma 11.2. For every x e M®, 



1=1 

The proof of the lemma is essentially the same. We then use the lemma to prove 
the uniqueness of the association scheme. 

Theorem 11.3. There is only one 4-class association scheme with the same size, 
valencies and intersection numbers as the association scheme of minimal vectors of 
the Eg lattice. 

The proof of the theorem involves, as before, the operator 



which turns out to be a projection to an 8-dimensional space. Again, the proof is 
essentially the same as for the Leech lattice. 

11.4. Inner product bounds. The inner product bounds use only the intersection 
numbers fLemma l7.4[ for which we needed the isomorphism of association schemes, 
deals with a case that does not occur in E'g). We get almost the same bounds as in 
Section [71 The proofs have to be slightly modified due to the fact that the minimal 
norm of the Es lattice is 2 instead of 4 for the Leech lattice. We get the following 
result: 

Lemma 11.4. Let u,v be nearly minimal vectors. 

(1) For {u, v)Ki2 we have 2 < (w, u) < 2 + (5/2)£. 

(2) For {u, -y) « 1 we have 1 - (5/2)e < {u,v) <l + (9/2)£. 

(3) For {u,v) fa we have ~7e < {u,v) < 7e. 




(3.48 • 10"-*) ^ + 30(7 + 240(1 + 2a)a < 4.4 • 10 



,-3 



N 
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11.5. A basis of nearly minimal vectors. The proof that there is a basis of 
nearly minimal vectors is completely analogous. Consider the basis of minimal 
vectors 



1 
1 






Vl/2 



1 

-1 
1 





1/2 






-1 
1 





1/2 







-1 
1 




1/2 








-1 
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1/2 









-1 
1 

1/2 










-1 

1/2 



\ 








1/2/ 



of £'8, where the row vectors are the basis elements. Then the corresponding vectors 
in A (via an isomorphism of association schemes) can be shown to be a basis for A. 



11.6. Computing a. One first proves an analogue of Lemma FlO.ll bv considering 
the vectors Wi = (0, . . . , 1, 1, . . . , 0) and Vi — (0, . . . , 1, — 1, . . . , 0) where the i and 
i + 1 coordinates are nonzero, for i € {1, 3, 5, 7}. 
The analogue of Lemma [10.21 is immediate: 

Lemma 11.5. Let vi,...,vs be orthogonal minimal vectors in the Eg, lattice. If 



0, then 



We then proceed to use these lemmas as before to prove bounds on a. 

Lemma 11.6. With the same notation for (3 and t as before, we have the following 
bounds: 

(1) If 13 = 2 and t < then we can take a — 1. 

(2) If (3 = 2 and t > then we can take a = 1/7. 

(3) If (3 — \ and t < then we can take a = 2/9. 

(4) If [3 — \ and t > then we can take a — 2/15. 

(5) If (3 — Q then we can take a = 1/20. 

In fact, the greatest possible value of a is 1/7, and that can be rigorously proved 
by computer calculations. (The linear programs are small enough that one can solve 
them using exact rational arithmetic. The Maple file E8seventh.txt contains the 
calculations.) However, the weaker bound of 1/20 will suffice. 



11.7. Local optimality of Eg. The proofs of perfection and eutaxy of Eg closely 
parallel those for the Leech lattice. The Gram matrix for the basis of Eg that we 
chose above is small enough that we can compute all its minors quickly. We find 
that for < p < 10"^, 

Dp>l- p2(7936 + 21162%/3p + 84256/9^ + 47300%/5p^ 
+ 74088p'' + 10290V7p-^ + 4096p^) 
> 1 - 7973p2. 
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We conclude from the above calculations that for every perturbation by ptij where 
max{|iij |} < 1 and J2i j — 0, we have some k such that 

^^'.J:^''^"^ < (1 - P/40)(l - 7973p^)-V8. 

This bound is strictly less than 1 when < p < 2.5 • 10~^. 

As before, we have one final issue to deal with: our perturbation may not satisfy 
J2i j ^ij'^ij — 0- We normalize as before by setting 

, _ S j j + rii j ^ 
'''^ l + A/8 

(the notation is as in the Leech lattice case), and need to check that \tij \ < 2.5- 10~^. 
From Lemma [11.41 we have \riij\ < 7e. Now the sum of the absolute values of the 
entries of S is 620, so |A| < 4340£. The maximum value of \sij\ is 2, so putting 
everything together we have 



77ij + Si J A/8 



^ 7. + 2.4340e/8 
- 1 - 4340e/8 



l + A/8 

This proves that A must be the same as Eg since it lies within the range of local 
optimality of Eg. 

This completes our new proof of the optimality of Eg: 

Theorem 11.7 (Blichfeldt, Vetcinkin). The Eg, root lattice is the unique densest 
lattice in M*, up to scaling and isometrics o/R^. 
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Appendix A. Computer calculations 

The computer files for checking our calculations are available from the arXiv . org 
e-print archive. This paper is available as math.MG/0403263. To access the auxil- 
iary files, download the source files for the paper. That will produce not only the 
files for the paper but also the computer algebra code. 

By far the most extensive use of computer calculations in this paper occurs in 
Subsection 14.11 The calculations are carried out in the file verifyf.txt, which 
consists of PARI code. PARI is a free computer algebra system designed for rapid 
number-theoretic calculations. See http://pari.math.u-bordeaux.fr for more 
information on PARI or to download a copy. Our PARI files all contain comments 
that should help make them understandable to those unfamiliar with PARI. 

First, we will explain how one proves the properties of the function / used in 
Subsection 14. 1[ then we will explain how it was constructed. Finally, we briefly 
discuss the verification of the other calculations in this paper. 

There are two auxiliary files for dealing with /: fcoeffs.txt contains coeffi- 
cients Co, ... , C803 which we will use to construct /, and roots . txt contains values 
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''Oj ■ • ■ , '"200 that are nearly roots (meaning the polynomial is very near at those 
locations, although in fact no real roots are nearby, except near ro). 
The file verifyf.txt carries out the following verifications. Let 

803 

(A.l) Mz) = J2c,z\Ll\z). 

1=0 

Define /: R^" R by f{x) = /o(27r|xp)e-'^l^l VlO^°°°. (The denominator of 
2^03000 jj^g^jjgg coefficients Ci integers, which is convenient.) The value of r used 
in Subsection 14.11 satisfies 2nr'^ = rg. 

We first need to check that f(x) < for |x| > r, which is equivalent to fo{z) < 
for z > rQ. In principle one could check this straightforwardly using Sturm's 
theorem, but that takes a tremendous amount of time for such a huge polynomial. 
Instead, we will use Descartes' rule of signs in a somewhat complicated way that 
may not appear a priori superior, but works overwhelmingly better in practice: the 
number of roots of a polynomial p in the interval (a, b) is at most the number of 
sign changes in the coefficients of 

and is congruent to it modulo 2. For (a, oo) one can simply use p{a + z). This 
result is sometimes known as Jacobi's rule of signs (see Corollary 10.1.13 in |RS[ 
p. 320]). 

We check using Jacobi's rule of signs that /o has no roots in ([7'2ool,oo) or 
{\ri], [ri+i\) for 1 < I < 199 (and also check that it does not vanish at the end- 
points), and that it has exactly one root in (0, [rij). Then we must check that it 
does not vanish on [[r^J , [r^]] with 1 < i < 200, and that its one root is less than 
ro (for that we simply compute the sign of /o(ro)). 

Dealing with the intervals [[r^J, [r^]] is a little more difficult. We use Jacobi's 
rule of signs to check that /q' has no roots on these intervals, so /q is monotonic 
and has at most one root in each. We then check that foivi — 10"^'''") > and 
/o(''i + 10^^'''") < 0, so the maximum of /o must occur within lO^'^^'^ of r^. However, 

|/^(r,± 10-350)1 < 103285^ 

from which it follows (by the mean value theorem and the monotonicity of /q) that 
fo{n)~ max fo{x) < 10^350 . ^03285 ^ ^02935 



Because |/o(ri)| < —10^^^, fo has no roots in [[r^J, [' 

To make the proof more efficient, we arrange the calculations to ensure that only 
integer arithmetic is used, so that PARI does not spend time reducing fractions to 
lowest terms. (That explains why the coefficients in (jA.l[) are multiplied by 

Dealing with / is similar, but of course it uses the polynomial 

803 

ho{z) =Y,{-lYc,tlLl\z) 

i=0 

instead of /o, and in this case plays the same role as ri, . . . , r2oo do (rather than 
being treated differently, as above). 
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All that remains is to check /(O) — /(O) = 1, which is true because /o(0) = 
/io(0) = lO"^"*^". These calculations suffice for the proof that no sphere packing can 
be more than a factor of 1 + 1.65 • 10~^° times denser than A24. However, Proposi- 
tion [4?3] requires more detailed information on the values of /. In verifyf .txt we 
check these inequalities using 

i=0 1=0 

for < z < 60, as well as rational upper and lower bounds for tt, to avoid having 
to deal with irrational numbers. 

Unfortunately the file verifyf . txt cannot address where / comes from (it takes 
far longer to locate / than to verify its properties). The construction is based on 
the numerical technique described in Section 7 of [CE| . which describes functions 
via forced double root locations that are then repeatedly perturbed until they reach 
a local optimum. Instead of this straightforward optimization algorithm we imple- 
mented a high-dimensional version of Newton's method (which locates a root of the 
derivative of /(2) as an implicit function of the forced double roots — this is valuable 
because /(2) is roughly proportional to |r — 2|). We arrived at 200 forced double 
roots, whose locations are specified in roots.txt. The next step of the method 
from |CE| also caused trouble: solving for the polynomial with those forced roots. 
Exact rational arithmetic was immensely time-consuming and produced huge de- 
nominators. Instead, we carried out high-precision floating point arithmetic and 
rounded the result to within 10^^°"°. That could ruin the sign conditions on / 
and /, by turning a double root into a pair of nearby single roots, so instead of 
solving for actual double roots we required that / and / should stay slightly on the 
correct side of 0. That does not greatly change the resulting bounds, and leads to 
the function / used in this paper. 

Lemma 14.41 requires the second largest amount of computation to check, al- 
though far less than Subsection 14.11 The file verifyg.txt and the auxiliary file 
gcoeffs.txt deal with it in a fairly straightforward way. 

The calculations required for (|9.2p and (|9.3p are dealt with in verifygram.txt. 

Finally, the file verifyrest.txt verifies all remaining calculations in the paper 
for the Leech lattice case. Many of the calculations in this file could be checked by 
hand, but that would be unpleasant and error-prone. Instead of PARI code, this 
file contains Maple code. Maple is less efficient but more fiexible, and it seemed 
easiest to use it for these calculations. 

Our computer calculations are completely rigorous, in the sense that they are 
carried out using exact arithmetic and thus avoid traps such as round-off error. 
Nevertheless, we cannot eliminate the possibility of a hardware error or a bug in 
software beyond our control, such as the operating system or computer algebra 
program. Under these circumstances, even assuming our programs are correctly 
written they may return incorrect results. This possibility appears quite unlikely, 
and we do not consider it a serious worry. However, we have addressed it by 
asking Dimitar Jetchev to translate our programs into the Magma computer algebra 
system. His translation is contained in the file magmacode.txt. Using it, we have 
independently checked our calculations on a different type of processor, a different 
computer algebra system, and a different operating system. 
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We have also documented our calculations for the Es proof. The details can be 
found in the files ESverifyf .txt and ESrest .txt. We include fewer comments in 
these files than in the others, because their structure is parallel to the Leech lattice 
case. The file E8seventh.txt contains a proof that a = 1/7 works in the Eg, case, 
although we do not require that for the proof of Theorem 111.71 

Appendix B. Background 

In this appendix we collect brief definitions and descriptions of a few of the 
principal objects and techniques used in this paper. 

An even unimodular lattice A C M" is a lattice such that |A| — 1, {x^y) G Z 
for all x,y € A, and (x, x) € 2Z for all x € A. Such lattices exist only when n is 
a multiple of 8. Up to isometrics of R", the unique example when n = 8 is Eg, 
and there are 24 examples in R^'', among which the Leech lattice is the unique one 
containing no vectors of length \/2. See |CS2( p. 48] for more information. 

A spherical code of minimal angle (p in the unit sphere S"~^ is a collection 
C C S*"^^ of points such that {x,y) < cosyj for all x,y £ C with x ^ y. In other 
words, no two distinct points of C form an angle smaller than Lp centered at the 
origin. Spherical codes are to S*"^^ as binary error-correcting codes are to {0, 1}", 
or as sphere packings are to R". 

The most important approach to bounding the size of spherical codes is linear 
programming hounds (due to Delsarte [D] ; see Chapter 9 of [CS'2] for an exposition) . 
These bounds rely on the following property of the ultraspherical polynomials, 
which follows from Lemma [Ol ifA = n/2— 1, then for each finite subset C C S*""^, 

^ C^{{x,y))>Q. 

x,yGC 

Suppose fo, . . . , fd > 0, and that f{z) = X^iLo fi^ii^) satisfies f{z) < for 2; e 
[— l,cos(/3]. Then every spherical code C in S'"~^ with minimal angle ip has size 
bounded by \C\ < /(l)//o (assuming /o 7^ 0). The proof is simple: 

|c|/(i)-E/((^'^)) 

xec 

x.yeC 

d 

= E /o + E/' E c^a^^y)) 
> E/o 

X.yeC 

= |C|Vo. 

The most dramatic application of the linear programming bounds for spherical 
codes is the solution of the kissing problem in R* and R'^'*: how many unit balls can 
be placed tangent to a central unit ball so that their interiors do not overlap? This 
condition amounts to saying that the points of tangency form a spherical code with 
minimal angle 7r/3. In R^ the answer is 240, and in R^* the answer is 196560. The 
codes are formed from the minimal vectors in Eg and the Leech lattice, respectively 
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(of course the radius of the sphere involved differs in the two examples). Optimality 
follows from the linear programming bounds by taking 

in and 

in M.^'^. This most remarkable fact was discovered independently by Levenshtein 
[Lev] and by Odlyzko and Sloane [OSj . 

A spherical t-design in S"'~^ is a non-empty finite subset D of S"^^ such that 
for every polynomial / : M" ^ R of total degree at most t, 

In other words, the average of / over V equals its average over the entire sphere. 
The minimal vectors of the Leech lattice form a spherical 11-design, and those of 
Es form a spherical 7-design. 

A k- class association scheme is a set A together with a partition 

A"^ ^ AaU---UAk 

such that Ao = : x G .4}, {x,y) £ Ai iff {y,x) G Ai, and the following 

property holds. Fix £, m, and n in Z fl [0, k]; then for all x,y G A with {x, y) G Ai, 
there are the same number Pi{m, n) ol z £ A such that (a;, z) G Am and (j/, z) G An- 
(That is, Pi{m, n) depends only on £, to, and n, and not on x and y.) These numbers 
are called the intersection numbers of the association scheme. When £ = and 
m = n they are also called valencies. In the body of the paper we modify this 
notation slightly for the association scheme C\ (in a trivial way): we label the 
classes with the corresponding inner products in C24, and use these labels in the 
notation for the intersection numbers. 
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